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Flavor asymmetry is investigated in polarized light-antiquark distributions by a meson-cloud 
model. In particular, p meson contributions to Au — Ad are calculated. We point out that the 
g2 part of p contributes to the structure function gi of the proton in addition to the ordinary 
longitudinally polarized distributions in p. This kind of contribution becomes important at 
medium x (> 0.2) with small (~1 GeV^). Including N pN and N ^ pA splitting 
processes, we obtain the polarized p effects on the light-antiquark flavor asymmetry in the 
proton. The results show Ad excess over Au, which is very different from some theoretical 
predictions. Our model could be tested by experiments in the near future. 
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I. INTRODUCTION 

Light antiquark distributions are expected to be al- 
most flavor symmetric according to perturbative quan- 
tum chromodynamics (QCD). The next-to-leading-order 
effects contribute to the difference between u and d; how- 
ever, the contribution is tiny as long as they are esti- 
mated in the perturbative QCD region. Therefore, it was 
rather surprising to find the antiquark flavor asymmetry 
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u/d in Gottfried-sum-rule violation by the New Muon 
Collaboration (NMC) and in succeeding Drell-Yan 
and semi- inclusive measurements of the CERN-NA51 [|], 
Fermilab-E866/NuSea §, and HERMES @ collabora- 
tions . In particular, the E866/NuSea experimental re- 
sults played an important role in establishing the fla- 
vor asymmetry by clarifying the x dependence of u/d. 
This new experimental finding was a good opportunity 
of investigating a mysterious nonperturbative aspect of 
hadron structure. 

Various models have been proposed to explain the un- 
polarized flavor asymmetry. So far, meson-cloud type 
models are successful in explaining the experimental re- 
sults. For the explanation of these models and other 
ideas, the authors suggest reading the summary papers 
in Ref Since most theoretical papers are written 

after the NMC finding, the actual test of the proposed 
models should be done by predicting unobserved quan- 
tities. In this sense, the polarized light-antiquark flavor 
asymmetry should be an appropriate one. In fact, there 
are already several papers on this topic by phenomeno- 
logical hadron models in Refs. |6| ^, ^. It is particu- 
larly interesting to find that a meson-cloud model and 
a chiral-soliton model predict totally different contribu- 
tions to Am — Ad, although their results are similar in 
the unpolarized distribution u — d. 

The situation of the polarized antiquark distributions 
is not as good as the unpolarized one in the sense that the 
polarized whole sea-quark distribution itself is not well 
determined at this stage. Most parametrizations assume 
flavor symmetric polarized antiquark distributions, which 
are then determined mainly by the gi measurements. As 
a result, there is much uncertainty in the antiquark dis- 
tributions at small and large x |g, |l^. Although there 
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are some semi-inclusive data which could be sensi- 
tive to the light antiquark flavor asymmetry, they are 
not accurate enough to provide strong constraint for the 
polarized antiquark flavor asymmetry [^o[ H owever, 
the Relativistic Heavy Ion Collider (RHIC) |12| and the 
Common Muon and Proton Apparatus for Structure and 
Spectroscopy (COMPASS) jl^ experiments should clar- 
ify the details of the polarized antiquark distributions in 
a few years. It is the right time to investigate the anti- 
quark flavor asymmetry Au/Ad by possible theoretical 
models and to summarize various predictions. 

In this paper, we intend to shed light on the virtual 
meson model which has been successful in the unpolar- 
ized studies . The purpose of this paper is to extend 
the studies of the virtual p-meson contributions by Fries 
and Schafer in Ref. 0|. In particular, we point out that 
the g2 part of the polarized p contributes to the polarized 
flavor asymmetry in addition to the ordinary longitudi- 
nal part, which was calculated in Ref. 0. Because we 
show new g2 terms in this paper and because the situa- 
tion is still confusing in the sense that another p-meson 
paper |^ claims major differences from Ref. |Q in sup- 
posedly the same p contributions, the detailed formalism 
is shown in the following sections. The meson model was 
extended recently to a different direction in Ref. |^ by 
including tt — p interference terms; however, this paper 
is intended to investigate a different kinematical aspect 
within the meson model. 

The paper consists of the following. The formalism 
of p contributions to Am — Ad is presented in Sec. p . 
Meson momentum distributions are obtained in Sec. [H, 
and numerical results of Au — Ad are shown in Sec. IV. 
Our studies are summarized in Sec. ^ 



II. VECTOR-MESON CONTRIBUTIONS 

The cross section of polarized electron-nucleon scat- 
tering is generally written in terms of lepton and hadron 
tensors: 



da 



dELdVLL 



\Pe\ 



(q2)2 



Lf^^iPe, Se, q) W^^{pN, SN,q), 



(2.1) 



where a is the fine structure constant, E'^ and fig are 
the scattered electron energy and solid angle, and Pe, 
Pg, pn, and q are initial electron, final electron, nucleon, 
and virtual photon momenta, respectively. The electron 
and nucleon spins are expressed by Sg and sn with the 
normalization Sg ~ = —1. Throughout this paper, 
the convention —goo — 9ii — 922 = 333 = +1 is used 



so as to have, for example, pj^ 



Pn 



I'M- 



Furthermore, the Dirac spinor is normalized as u^u = 
Ee/m^ or En/ TUN, where E^ and Epj are electron and 
nucleon energies, and mg and rriN are their masses. The 



polarized lepton and hadron tensors are given by 

L^''(pg, 5g, q) = 2 [pj^p^"- + p:p'/ - (pg . - ml) g'^- 
-ie^^'^p'^m^qps,^], (2.2) 

Wf,^{pN,SN,q) = ^^(27r)^(5''(pAr + g-p^) 

X 

X <pN,SN\J^.mX >< X\MO)\pN,SN >, (2.3) 

where the factor et^'^P" is the antisymmetric tensor with 
the convention £9123 = +1- 




FIG. 1: Virtual vector-meson contribution. 

Next, we consider the process in Fig. ^, where the 
nucleon splits into a virtual vector meson and a baryon, 
then the virtual photon interacts with the polarized me- 
son. Because scalar mesons do not contribute directly to 
the polarized structure functions due to spinless nature, 
the lightest vector meson, namely p, is taken into ac- 
count in this paper. In future, we may extend the present 
studies by including heavier vector mesons. As the final 
state baryon, the nucleon and A are considered. Express- 
ing the VNB vertex multiplied by the meson propagator 
as JvNB{k, sv ,Pn , sn ,Pb, sb) and calculating the cross 
section due to the process in Fig. |^, we obtain 

dE'jn', ~ |pg| (g2)2^ ^P-^'-^^> 

X \JvNB?<Ksy\J,mX> 



(27r)3 2tiEb 



X < X\J^{Q)\k,sv > {2TTY5^{k + q-p^) 



(2.4) 



Here, k and sy indicate the meson mom entu m and spin. 
This equation has the same form as Eq. ( |2.lD . Therefore, 
the last part is identified with a vector-meson contribu- 
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tion to the nucleoli tensor: 

, . I (PpB 2mvmB 



(27r)3 Eb 
\JvNB\^W^]:Hk,sv,q), (2.5) 



where my is the meson mass, and the meson tensor is 
defined by 

X < k,sv\Jy,{fi)\X> < X\J^{0)\k,sv > ■ (2.6) 

In this way, the vector-meson contribution to the nu- 
cleon tensor is expressed in terms of the VNB vertex 
and the meson tensor. Because we are interested in me- 
son effects on the polarized parton distributions in the 
nucleon, we try to project the gi part out from the nu- 
cleon tensor. The definition of the gi and g2 structure 
functions is given in the asymmetric part of the nucleon 
tensor: 



91 



PN ■ q 



+ [PN ■ qsN - SN ■ qpn) 
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[pN ■ qf 



(2.7) 



In order to discuss each structure function separately, a 
projection operator 



2pN ■ q 
is then applied to give 



l£ qa SNf3, 



(2.8) 



P"-" W^^ipN, SN, q) = j—^ [q' + {s^ ■ qf] g, - 7^ g,. 



Here, 7 is defined by 



7 



{pN ■ q) 



2 4x^m^ 



(2.9) 



(2.10) 



with = —q^. In the same way, gi and g2 structure 
functions of the vector meson are defined in the asym- 
metric part of the tensor. Operating the projection also 
on the meson tensor, we obtain 

P'^'Wl^V^'ik, sv,q)^'^[A, gYik, q) + g^{k, q)], 

(2.11) 

where Ai and A2 are given by 



Ai 



PN ■ qk ■ q 



(sN ■ qsv ■ q - q sn ■ sy), 



(2.12) 



mNmvq^ , , , 

A2 = 7-. [SN ■ ksv ■ q - k ■ qsN ■ sv). (2.13) 

PN ■ q[k ■ qy 



From Eqs. (2^), (2^), and ( 2.11 ), the meson contri- 
bution to the nucleon structure functions becomes 



'TV 



(pN ■ qy 



[q^ + (sN ■ qf]giiPN,q) - i^g2{pN,q) 



\Jvnb\^ 



cPpB 2mNmB \ ^ 
X [A,gYik,q)+A2g^{k,q)-_ 



(2.14) 



Then, the above integration variables (p^, p^, pg) are 
changed for the meson momentum fraction y, the trans- 
verse momentum k± of the meson, and the angle (j) be- 
tween k±_ and the transverse spin vector of the nucleon 



y 



k ■ q 

PN ■ q 



fc • Siy = kj_TN cos ( 



(2.15) 



with T]\f = |s^|. Then, the meson contribution is ex- 
pressed as 



i-N 



{pN ■ q)' 



[q^ + [sN ■ qy]giix,Q^) - g2ix,Q^) 



' ^ [B,{y) gYix/y, Q^) + B2{y) gX{x/y, Q^)]. 

X y 

(2.16) 



The upper limit of the y-integration range is taken as 1 by 
considering the vector-meson mass smaller than the nu- 
cleon mass. However, one should be careful in extending 
the present studies to other mesons with larger masses. 
The meson momentum distributions are expressed as 



BiAy) 



dk: 



2Tr 



\pN\mN TUB dy' 
{2^fEB dy 



y 



"'0 
X \JVNB?A^^2 , (2.17) 

where y' is the longitudinal momentum fraction defined 
in the meson momentum 



k^k±+ y'pN- 



(2.18) 



In the infinite momentum frame \pn\ ^ 00, y and y' are 
related by 



y 




y^m% 



(2.19) 



Because time-ordered perturbation theory is used for the 
reaction in Fig. || as explained in Sec. |l|, the vector 
meson is taken as an on-shell particle: my in the 

above derivation. The partial derivative dy' /dy can be 
calculated from this expression. In the infinite momen- 
tum frame, the momentum fraction y' has to satisfy the 
kinematical condition < y' < 1, namely the meson V 
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and the baryon B should move in the forward direction. 
The maximum transverse momentum is given by 



'AT 



r 



( + + 1) (a/1 + 72 + 1 - 2j/) - m 



(2.20) 

Practically, it does not matter to take the upper bound 



{k'^)max — > oo in Eq. (2.17) at small-a; where the anti- 



quark distributions play a major role, because {k'j_)max 
(~ Q^{l — y)/x^ ^ tt^'n) is beyond the vertex momentum 
cutoff region discussed in Sec. III. The contribution to 



the integral between {kj_)max and fc^ = oo is extremely 
small in general. Furthermore, the upper bound becomes 
ik]_)max — > oo in the limit oo, and it is consistent 

with the previous publications [0, |^. In this way, the 
meson contribution is expressed in terms of the meson 
structure functions convoluted with the meson momen- 
tum distributions in the nucleon. 

Using the integration variables y, k'^ , and 0, we express 
the coefficients Ai and A2 as 



yy' m^N 



-(VT 



1) 



7^TArAy COS (j) 



ki 



A. = 



9 9 

y'^m% 



yrriN 

— AjvAy + Tat AyCOS ( 



(VT+7^-1) 

(2.21) 



in the limit \pn\ 00. Detailed calculations indicate 
that (f> dependence can be extracted out from another 
part of the integrand in Eq.(2.17) as 



As 



Tat cos 4> Cj^' . 

(2.22) 



Then, after the (j) integration, Eq. ( 2.17 ) becomes 



Bi{y) 



Ay [\Nl^l{y)-rlf^-{y)], 



Mr 

B2{y) = V [-Ajv f^l (y) + tI f^- (y)]. 



Av 



(2.23) 
(2.24) 



where Av dependence is explicitly denoted in meson mo- 
mentum distributions, which are defined by 



fi^y) 



fi^'iy) 



f2l{y) 



dklj^nCj 



yym% 



yrriN 



dkhHnCt^^, 



(2.25) 
(2.26) 

(2.27) 



f2^{y) 



k±'my 
y'y^m% 



(v/T+^-1). 

(2.28) 



Because the functions f^^, f^]^, and f^f are propor- 
tional to 7^, they van ish i n the limit Q'^ ^ 00. As it 
is obvious from Eq. ( ^.H ), it is necessary to consider 
both longitudinal and transverse polarizations for the nu- 
cleon in order to extract the gi part. In addition, the 52 
structure function of the meson contributes. The func- 
tion J^^ (y) is the ordinary meson momentum distribu- 
tion with the momentum fraction y in the longitudinally 
polarized nucleon. The function /^-p (y) is the distribu- 
tion in the transversely polarized nucleon. On the other 
hand, (y) ^^'^ /2T (2/) ^'^^ distributions associated 



with g2 of the vector meson. Expressing Eq. (2.16) in 
terms of the nucleon and meson helicities, Aat and Ay, 
we obtain 

(A'^ - rjj 7') 51 {X, Q2) _ ^2 ^^(^^ g2) 

= J2xv /iL(y) ~ t'n fiTiy)}9Y{^/y, Q') 

\ J X y 

+ { - Aat /2L(y) + 4 f2T{y)} g^ix/y, Q^) ] . (2.29) 

Combining the longitudinal polarization Aat = 1 {tn — 
0) with the transverse polarization tn — I (Aat = 0), we 
can extract the gi part as 



gi{x,Q 



'dy 

y 



1+7^ 

X [{AAi(y) + A/iT(y)}.gr(x/y,Q2) 
- {Af2Liy) + Af2T{y)} gYix/y,Q'y 



(2.30) 



where the functions /S.fY^ {y) with i=lL, 2L, IT, and 
2T are defined by 



A/.(y) = /,"^'=+^(y)-/, 



fAv = -l 



(2.31) 



The g2 part is obtained in the same way as 
1 r'dy 



92{x,Q ) = 



1+7" 



y 



X [{-A/ii(y) + A/iT(y)/7'}5r(^/y, 
+ {AM(y) - A/2T(y)/7'} 9X(x/v,Q^)] 



(2.32) 



In the limit — > 00, namely 72 ^ 0, only the momen- 
tum distribution A/ii(y) remains finite, and Eq. ( |23o| ) 
agrees with t he ex pression in Ref . 0| . 

In Eq. (2.30), there are additional terms associ- 



ated with g2 of the meson. For discussing these 52 
type contributions to Au — Ad, g2 is approximated by 
the Wandzura-Wilczek (WW) relation ||l^ by neglecting 
higher-twist terms: 



^gX{x,Q^ 



—91 (2/, Q ) 



(2.33) 
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Then, providing the leading-order expression for gY , we 
have 



92 



+AqT^\.,Q^)]. (2.34) 
The above WW distributions are defined by 



J X 



(2.35) 



and the same equation for AqY^^^\x, Q^). From these 
equations, we obtain a vector meson contribution to the 
polarized antiquark distribution Aq^ in the proton as 

AC^^(x,Q^)=-^ 

^ + 1 Jx V 
X [{A/ii(y) + A/iT(y)}AC(x/2/,Q2) 

- {Ah^iy) + Ahriv)] Aqf^'^^) (x/y, Q') ] . 

(2.36) 

If this kind of vector-meson contribution is the only 
source for the polarized flavor asymmetry, the Au — Ad 
distribution is then calculated by taking the difference 
^qY=u^ - Ag]Lrf^ in the above equation. 



III. MESON MOMENTUM DISTRIBUTIONS 

In order to estimate the meson contributions numeri- 
cally, it is necessary to calculate the momentum distribu- 
tions AfY^ (y) of the meson. We calculate them by con- 
sidering the vector-meson creation processes TV — > VN' 
and N VA through the interactions 



VvNN = 4>v ■ T FvNN{k) u{p^, s'^) e^* 
fv 



ZruN 



u{pn,sn), (3.1) 



VvNA = 4>v ■ T FvNA{k) UuipA, sa) ^^J^^ Jalfi 



my 



(3.2) 



where u{p]\j,sn) is the Dirac spinor, ?7^(pa,sa) is the 
Rarita-Schwinger spinor, and e'^ is the polarization vec- 
tor of the vector meson. The VNN and VNA coupling 
constants are denoted as gv, /y, and fvNA, and form 
factors are denoted as FvNNik) and Fyj^Aik)- Isospin 
dependence is taken into account by the factor (j)y ■ T, 
and it is defined in terms of a reduced matrix element 
and a Clebsch-Gordan coefficient lll6|, ITtII 



<B\<j)*yT\N> = {-l) 



<Tb\\T II Tn> 



V2Tb+T 
X < TjvMjv : 1 - Mv\TbMb > (3.3) 



with < 1/2 II f II 1/2 >= VQ and < 3/2 || f || 1/2 >= 2. 
Here, Tpf and Tg denote isospins of the nucleon and the 
baryon, respectively, and Mj^ and AIb are their third 
components. 

From these vertices, the meson momentum distribu- 
tion fniy) can be calculated together with the baryon 
distribution fB{y)- They are supposed to satisfy the rela- 
tion fniy) — fsi^ — y) because of charge and momentum 
conservations. However, it is known that the covariant 
calculation could violate this relation because a deriva- 
tive coupling introduces off-shell dependence. A possible 
solution 1^, |l^ is to use the time-ordered perturbation 
theory (TOPT), instead of the covariant formalism. Al- 
though the four-momentum conservation is satisfied at 
the VNB vertex in the covariant formalism, the energy 
is not conserved in the TOPT |2^ . If there is no off-sheU 
dependence at a vertex, the TOPT agrees certainly with 
the ordinary covariant theory by collecting all the time- 
ordered diagrams. However, the off-shell dependence due 
to the derivative coupling complicates the problem. It 
leads to a fre edom in d efining the vertex momentum K 
in Eqs. (3.1) and (3.2). The following two possibilities 
are considered in Refs. |0, |l8|: 

(A) k = k = {Ev,k), 

(B) k^k=pN-pB = {EN-EB,k), (3.4) 



my + k'^. There is another off-shell de- 



where Ey - 

pendence from the vertex form factors, and it is discussed 
in Sec. ITv. 



From the VNN interaction vertex in Eq. (3.1), we 
obtain 



9l{-k" 



+ 2 (pat • ep'j^ ■ e* +PN ■ £*p'n ■ e) + 2mNmvSN ■ sy} 
+ gv fv { — 2k ■ K + 2mNmySN ■ sy 

~ 2p'j^ ■ k SN ■ s' - k ■ ek ■ e* - k ■ e*k ■ e} 



+ -p^ {{k^+k-ek- e*){k^ - 4m?^) + ApN ■ kp'j^ ■ k 

- 2K'^{pn ■ epN ■ £* +PN ■ £*Pn ■ e) + ^m^K'^SN ■ s' 

— irriNmyp'j^ ■ Ksjy ■ Sy — 2m]ymyK'^SN ■ sv 

+ 2{pN ■ k +p'i^ ■ k){pN ■ sk ■ e* +PN ■ £*k ■ e) 



2pN ■ k{k ■ ek- e* + k ■ e*k ■ e)} 



where s"^ , Sy, Sy are defined by 



TUN 



^y = S^''"^ st Sa^kp, 

my 

3y = £^ ^ £y £a kf-j. 

my 



(3.5) 

(3.6) 
(3.7) 
(3.8) 
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The meson polarization vector is given by 



mv 



mv{Ev + mv) 



(3.9) 



where Xy is the meson hehcity. The spherical unit vector 
exv is defined in the frame with the z' axis parallel to k: 

e±i^T^ix' ±iy'), eo = z'. (3.10) 



In the same way, the VNA term is calculated from Eq. 
(Fl) as 



Z 0.06 

CL 

< 0.04 
0.02 



Aa 



VNA 



(fc) 



f? 



VNA 



6 rriNm^my 

- 2m\pN ■ KpA ■ K + 2mNm\{K^ + K ■ eK ■ e*) 

+ 2pN ■pA{rnAK'^ + m\k ■ eK ■ e* + K'^pa ■ epA ■ e* 
- PA ■ K{pA ■ek-e*+pA- e*K ■ e) - {pa ■ K)^} 

+ "'*a{ - Pn ■ K{pA ■ eK ■ e* + PA ■ e*K ■ e) 
- PA ■ K{pN ■ eK ■ e* +pn ■ £*K ■ e) 
+ K'^ipN ■ epA -e* +PN ■ £*PA ■ e) } 

+ niNSjy ■ s'{2m]y{pA ■ KY + m\{mN + mA)K'^] 

— mjsimvSM ■ sv{2{pA ■ K)"^ + m^if^} 

— AlTLNlTlAmvPA ■ K SN ■ sv 

+ 2m1fmAmv{K'^SN • Si — Pa ' Ksn ■ h 
where Si and st^ are defined by 



0.2 0.3 0.4 0.5 D.6 0.7 0.8 0.9 

y 

FIG. 2; Meson momentum distributions A/iL(y) and 
A/2l(2/) from tlie pNN process. The solid and dashed curves 
are obtained at — 1 and 2 GeV^, respectively, with x—0.2. 
The isospin factors are taken out from the distributions as 
explained in the text. 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



FIG. 3: Meson momentum distributions Afiriy) and 
A/2r(t/) from the pNN process. The solid and dashed curves 
, (3.11) are obtained at — 1 and 2 GeV^, respectively, with x=0.2. 



'1 ~ 2 



rriN ruA my 



{pN ■ eel- PN ■ £*£u)PNa kp , 

(3.12) 

J^""'' {PN -eel-pN- e*e^) pAa Kp . 

(3.13) 



Therefore, Jvnb is expressed as 
JvNB 



VvNB- 



(3.16) 



Using Eqs. f^J), (|3.1lD , and ( 3.16| ) together with Eqs. 
and (|2.25D-(P^, we can calculate the meson mo- 



In Sec. U, the term Jvnb is defined by the vertex 
VvNB multiplied by the meson propagator. The propa- 
gator is the addition of two time-ordered terms. However, 
only the first one remains finite in an infinite momentum 
frame pn — > oo- 



1 



1 



2Ey{En — Ey — Eb) 2Ey{EB — Ey — En 

1 



where m\rg is defined by 



"T-ys = {k+PB^ 



my + kj_ ^ + kj_ 

y' 1-2/' 



(3.14) 



(3.15) 



mcntum distributions. The actual calculations are par- 
tially done by a Maple program. Obtained expressions 
are rather lengthy, so that the results are written in Ap- 
pendix. 

The momentum distributions are numerically calcu- 
lated by using the expressions in Appendix. However, the 
derivation of these analytical expressions is complicated, 
and it could easily lead to a calculation mistake. In order 
to avoid this kind of failure, we calculated the momentum 
distribution s nu meric ally i n an independent way directly 
from Eqs. ( j3.5| ) and (3.11), and we confirmed that they 
indeed agree on the results in Appendix. 

We show the numerical results in Figs. |^-|^ for the 
vertex momentum (B), which is the preferred choice ac- 
cording to Ref. [Esj. In addition to the variable y, the 
distributions depend also on x and Q^. The distributions 
are calculated at a; = 0.2 with and 2 GeV^ for the 

solid and dashed curves, respectively. Because of the x 
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' -,11, 





0.2 0.3 0.4 0.5 



0.7 0.8 0.9 



FIG. 4: Meson momentum distributions A/i_l(i/) and 
A/2L (y) from the pA^ A process. The sohd and dashed curves 
are obtained at — 1 and 2 GeV'^, respectively, with x—0.2. 




0.2 0.3 0.4 0.5 



0.7 0.8 0.9 



FIG. 5: Meson momentum distributions Afixiy) and 
A/2t(j/) from the pNA process. The solid and dashed curves 
are obtained at Q^=l and 2 GeV^, respectively, with x—0.2. 



dependence, the unphysical region y < x is not shown in 
these figures. In Figs. |^ and || the distributions due to 
the pNN process are shown. In Figs. ^ and ||, the dis- 
tributions due to the pNA process are shown. In these 
distributions, the isospin factors are taken out from the 
distributions, so that the distributions Af{y)/\(j)y • 
are actually shown, and there are differences of factors 
3 (in pNN) and 2 (in pNA) from those in Ref. 0. 
The coupling constants are taken as 5^/(47r) = 0.84, 
fv = G.lgv, and f^^Ji^n) = 20.45 Q. 

In spite of the positive A/ii(y) and A/2L(y) distri- 
butions from the pNN process in Fig. ^ the distribu- 
tions from the pNA are mainly negative. In the un- 
polarized case, the meson momentum distributions are, 
of course, positive. However, because the distri butio n 
A/(y) is defined by the helicity difference in Eq. ( 2.31 ), 



it becomes either positive or negative depending on the 
helicity structure at the VNB vertex. For the meson 
angular momentum state (.z , the VNB vertex amplitude 
is proportional to fc^' and higher order terms. Since the 
momentum k± is in general much smaller than the nu- 
cleon mass, the vertex amplitudes with £z — contribute 
dominantly to A/(y). There is only one pNN amplitude 
with iz ^ and Ay ^ 0, and it has the helicity structure 
Ay = +1 and A^ = —1/2, where the initial helicity is 



fixed at Aat = +1/2. This fact indicates that fil'~^^{y) 
and f2L~^^iy) certainly larger than fiL~~^iy) and 
fiL~~^{y)j respectively, which results in the positive dis- 
tributions A/ii(y) and A/il(j/) from the pNN process. 
On the other hand, there are two amplitudes with £z = 
and Ay 7^ in the pNA process, and they have helicity 
states Ay = —1, Aa = -f3/2 and Ay = -fl, Aa = —1/2. 
Actually calculating these helicity amplitudes, we find 
that both amplitudes depend much on the momentum 
choice, namely (A) or (B), at the pNA vertex. Therefore, 
fiL~^^iy) ^^'^ f2L~^^(y) either larger or smaller 
than fiL^~^iy) and f^l^'^iy), respectively, depending 
on the momentum choice. In the prescription (B), the 
positive helicity distributions fil~'^^{y) and f2L~^^iy) 
are mostly smaller, so that Afi^iy) and A/2L(y) become 
negative distributions in the wide a: region. However, the 
situation is opposite in the prescription (A), where the 
distributions are mostly positive. 

As expected, the distributions A/ii(y) in Figs. |2|and^ 
are the dominant ones and they are almost independent 
of However, A/2L(y), A/it(?/), and A/2T(y) are 
roughly proportional to so that these new contri- 

butions become more important as becomes smaller. 
Figures ||— |] clearly show this tendency. The transverse 
distributions A/ij'(y) and A/2t(2/) are an order of mag- 
nitude smaller than the longitudinal ones A/ii(y) and 
A/2l(2/). Therefore, the major correction comes from the 
distribution Af2L{y), which is almost comparable mag- 
nitude with A/ii(?/) in Figs. ^ and ^. Because the cor- 
rection terms are proportional to 7^ = Amj^x'^ / Q'^ , they 
are small contributions in the kinematical range x < 0.05 
with >1 GeV^. However, their effects become more 
pronounced as x becomes larger. 



IV. RESULTS 

For calculating the Am — Ac? distribution numerically, 
we need the polarized antiquark distributions in p, the 
isospin factors, and the vertex form factors. The p- meson 
parton distributions arc not known, so that the same 
prescription is used as the one in Refs. [0,^. Considering 
a lattice QCD estimate [Q, the polarized valence-quark 
distribution is assumed as 



AVp{x,Q^)=0.6V^{x,Q^), 



(4.1) 



at Q^ = l GeV^. The distribution in the pion is taken from 
the GRS (Gliick, Reya, and Schienbein) parametrization 
in 1999 [g3|. The charge symmetry suggests the relation 
for the valence-quark distributions: 



(An) 



val 
P~ 



2{AuY/ = 2(AJ)7' = AV, 



val 

p» 



(4.2) 



For the sea-quark distributions, they are assumed to be 
flavor symmetric. Then, we obtain the Au — Ad distri- 
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butions in the p meson: 



(Au - Ad)pO = 0, 
(Au - AJ)p- = +Ayp. 



(4.3) 



For the 52 part of p, the Wandzura-Wilczek relation is 
used as discussed in Sec. 



J X y 



'1- 

(4.4) 



Both the valence-quark distribution and the WW distri- 
bution are shown at Q2=1 GeV^ in Fig. | 



> 




FIG. 6: Assumed polarized valence-quark distribution in the 
p meson and the WW distribution at — 1 GeV^. 



Necessary isospin factors are calculated from Eq. (3.2) 



as 



|<n|0;+.Tb>|2^2, 
<A"|0:+.r|p>|2 = i/3, 



|<A^ 



:--T\p>\'^i. 



(4.5) 



Using Eqs. ( |2.36D , (^), and (|4j), we obtain 

{ A/ii + A/iT }(S>{Au- Ad), 



(Aw - Ad)p^pB = J2 

p,B 



{ A/2L + A/2T } ® (Au - Ad)^ 



WW 



2A/f 



L + IT + 3 



pNA 
L + IT 



' AK 



2 A/: 



A/: 



2L+2T r g ^J2L+2T 



Av: 



WW 



(4.6) 



where (8) indicates the convolution integral in Eq. (2.36) 



a(E)b ■ 



1 



^ + T Jx y 



1 dy 



-^a{y)b{x/y). (4.7) 



The meson momentum distributions A/^^^j, and 
A/2L+2T ^'^^ defined by extracting the isospin factors: 



A/, 



pWB _ AfiL + AfiT 
iL+iT ~ 



2. 



(4.8) 



The expression of Eq. ( 4.6D may seem to be different 
from Refs. ^ even in the limit 00; however, it is 

just the matter of the definition of the meson momentum 
distributions. They included the isospin factor 



\<n\cl>l^-T\p>\^ + \<p\c^\-T\p>\ 



(4.9) 



in the distribution Af{y) for the pNN process and the 
factor 



I<A0 



■T\p>\ 



|<A^ 



T\p>\ 



|<A++|C_.Tb>|2 = 2 (4.10) 



for the pNA. Therefore, our expression certainly agrees 
on those in Refs. @, |l at 00. 

The remaining quantities are the vertex form factors. 
They are roughly known from the studies of one-boson- 
exchange potentials (OBEPs); however, a slight change 
of the cutoff parameter could result in a large difference of 
antiquark distributions. Furthermore, there is an issue of 
the charge and momentum conservations for the splitting 
process [g4| if a i (= {pn —PbY) dependent form factor is 
used. A possible solution is to use the t dependent form 
factor multiplied by a w dependent one p5| . For this 
purpose, it is more convenient to take an exponential 
form factor so as to become the additional form t A- u 
within the form factor: 



FpNN{k) = FpNA{k) = exp 



"-VB 



2A2 



(4.11) 



where my^ is defined in Eq. ( 3.15| ), and the cutoff 
parameter Ae is taken as Ae=l GeV in the following 
numerical results. In Ref. |]l8|, the cutoff parameters 
are obtained by fitting baryon-production cross sections 
pp BX: Af^^^l.lO GeV and Ap^^=0.98 GeV. How- 
ever, the parameters are not well determined in general. 
We discuss the dependence on this cutoff value at the 
end of this section. The form factors are the same as the 
ones in the previous publications iQ, ^, so that we could 
compare our results with theirs. 

Using these form factors and the parton distributions 
in p, we obtain the pNN and pNA process contributions 
to the Au — Ad in the nucleon. In Fig. ^ the IL, 2L, 
IT, and 2T type distributions from the pNN process are 
shown at (5^=1 GeV^ together with their total. The ordi- 
nary IL term is the dominant contribution; however, the 
2L term becomes important at a; > 0.3. It is as large as 
the IL distribution in the medium x region although it is 
fairly small at a: < 0.05. The IT and 2T distributions are 
very small in the whole x range. Because p p^n is the 
only contributing process in which the valence d distri- 
bution in p^ plays the main role, the pNN contributions 
are negative in the Au — Ac? in the nucleon. 
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V / total 


Q-= l.OGeV- 
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FIG. 7: Aii — Ad distributions from the pNN process at Q^ = l 
GeV^ The IL, 2L, IT, and 2T type contributions and their 
summation are shown. 
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FIG. 10: A-u - Ad distributions from the pNN and piVA 
processes at — 1 GeV^ for the vertex momentum choice 
(A). 
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FIG. 8: All — Ad distributions from the pN /S. process at Q =1 
GeV^. The IL, 2L, IT, and 2T type contributions and their 
summation are shown. 
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FIG. 9: Au - 

cesses at 



Ad distributions from the pNN and pNA pro- 
1 GeV^ 



Each term contribution has almost the same tendency 
in the pNA process as shown in Fig. the IL term is 
the major one and the 2L term provides some corrections 
depending on the x region. There are two contributing 
processes, p p+A*' and p — s- p^A~, and the isospin 
factor is three times larger in the latter one. This fact 
may seem to indicate that the pNA processes provide a 



positive contribution to Au — Ad in the nucleon due to 
the valence u distribution in . This kind of explanation 
is certainly valid in the unpolarized flavor asymmetry 
1^, However, this is not the case in Fig. ^ where 
the IL and IT distributions are mostly negative. This 
misleading result comes from the helicity structure at the 
A'' pA vertex. Although the helicity difference Af{y) 
is positive for the pNN, it is negative for the pNA in 
the case (B) as explained in Sec. HI. Therefore, the 
pNA contribution becomes also negative for the Au — Ad 
distribution. 

Next, the pNN and pNA contributions are compared 
in Fig. ||. The magnitude of the pNA contribution is 
very small compared with the pNN one in (B). From 
Figs. I and |, we find that the magnitude of AfP'^'^{y) 
is already three times smaller than A/''^^(y), and the 
pNA contribution is further suppressed by the isospin 
factor (2/3)/2=l/3. Therefore, the overall magnitude 
becomes much smaller. 

As discussed in Sec. Ill, we may have another vertex 
choice (A) instead of {B)7 In showing the numerical re- 
sults so far, the model (B) has been used. We show the 
choice (A) results in Fig. |lO|. It is obvious that the dis- 
tributions depend much on this vertex choice. There are 
two major differences from Fig. ^. One is that the order 
of magnitude is much smaller in the pNN distribution, 
and the other is that the pNA distribution becomes pos- 
itive. These are due to the difference of helicity structure 
at the vertices between (A) and (B). 

We also discuss the vertex cutoff dependence. The ver- 
tex cutoff has been taken as Ae=l GeV in this section; 
however, it is well known that calculated antiquark dis- 
tributions are very sensitive to the cutoff value j26j. In 
the present paper, the exponential form factor is used 
instead of dipole or monopole form factor, which is more 
popular in the studies of OBEPs. The cutoff parameters 
of different form factors could be related by 

Ai = 0.62A2 ^ 0.78%/2Ae, (4.12) 

where the monopole and dipole parameters are defined 
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FIG. 11: Cutoff dependence of Am — Ad is shown at Q^=l 
GeV^ by taking Ae=0.5, 1.0, and 1.3 GeV. 



cess over Ad and the order of magnitude of Au — Ad 
is large (~ 0.4) at a;=0.2. Although the sohton models 
and the meson-cloud models obtain very similar distri- 
butions for the unpolarized u — d, it is interesting to find 
opposite prediction for the polarized distributions. The 
physics reason for the difference is not clear at this stage. 
In any case, the distribution Au — Ad should be clari- 
fied experimentally in the near future by W production 
process @ and semi-inclusive experiments |l^ . Further- 
more, there is a possibility to use the polarized proton- 
deuteron Drell-Yan process in combination with the 
proton-proton reaction. Until the data will be taken, the 
theoretical predictions should be discussed in details for 
comparison. 




0.5 0.6 0.7 0.8 0.9 



1 1.1 1.2 1.3 1.4 1.5 



Ae (GeV) 

FIG. 12: Cutoff dependence of \Au — Ad\ is shown at a; = 0.2 
and Q^=l GeV^ as a function of the parameter Ae. 



by the form factors 



(1) 

VNB 



(fc) = 



1 - "'w/Al p(2) 

^1 



VNB 



2^2 • 



Mi) 

(4.13) 



Even in the well investigated pion-nucleon coupling, 
the monopole cutoff parameter Ai ranges from 0.6 GeV 
to 1.4 GeV in quark models and OBEPs It roughly 
corresponds to Ae from 0.5 GeV to 1.3 GeV according 
to Eq. (4.12). We show the Au — Ad distributions for 
various cutoff parameters, Ae=0.5, 1.0, and 1.3 GeV in 
Fig. |l^ for the prescription (B). In the unpolarized dis- 
tribution u — din Ref . |^ , it is fortunate that the cutoff 
dependence is rather small because of the cancellation 
between the ttNN and ttA^A processes. However, the 
pNN is the dominant contribution, as shown in Fig. H, 
in the present polarized studies of (B), so that the overall 
magnitude is much dependent on the cutoff parameter. 
There are orders of magnitude differences between the 
three curves. Next, fixing x at 0.2, we show the cutoff 
dependence in Fig. |l2|. In fact, there is four orders of 
magnitude variation from Ae=0.5 to 1.5 GeV. Therefore, 
an accurate determination of the cutoff parameters is a 
key for a reliable meson-cloud prediction. 

There have been studies in chiral soliton models 
They predict very different distributions, namely Au ex- 



V. CONCLUSIONS 

The p meson contributions to the polarized antiquark 
distribution Au — Ad have been investigated. In particu- 
lar, we pointed out that the (72 part of p contributes to the 
polarized distributions in the nucleon. We obtained the 
extra contributions denoted as 2L, IT, and 2T in addi- 
tion to the ordinary one {IL). Although the extra terms 
are small in the small x region (x < 0.05), the magnitude 
of the 2L term becomes comparable to the ordinary one 
in the x region x > 0.2. The contributions are impor- 
tant in the kinematical region of medium x with small 
Q^. The obtained Am — Ad is very sensitive to the cutoff 
parameter. The model should be investigated further in 
order to compare with future experimental data. 
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APPENDIX A: ANALYTICAL EXPRESSIONS 
OF MESON MOMENTUM DISTRIBUTIONS 

In the limit — )■ 00, the following meson momentum 
distributions agree on those by Fries and Schafer (FS) 
[Q with a minor misprint in a pN A term. The situation 
of the momentum distributions is somewhat confusing in 
the sense that Cao and Signal (CS) [^| pointed out two 
major differences from Ref. |^ although the formalism is 
exactly the same except for interference terms. Accord- 
ing to Ref. all the gpNNfpNN {gvfv in our notation) 
terms should be replaced by —gpNNfpNN, and the mo- 
mentum (B) results for pNA agree on those of (A) in 
Ref. instead of (B). 
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In spite of their claim, we believe that the FS results 
are right with the following reasons. We also checked 
the helicity amplitudes in Ref. [|8|, which is referred 
to as Jiilich in the following. In addition to obvious ty- 
pos, our results differ from the Jiilich expressions. First, 
complex conjugate should be taken if their expressions 
are given for the process N — > pB as indicated in their 
appendix. Second, /„ terms have different sign. If the 
Jiilich amplitudes were written for Np B oi B ^ pN, 
we would agree on their expressions. Depending on the p 
momentum direction, the /„ term in Eq. ( p.l| ) becomes 
either positive or negative, which could lead to the dif- 
ferent sign of the gvfv terms. However, it is obvious that 
the outgoing meson is considered in the formalism. Fur- 
thermore, taking summations of the helicity amplitudes, 
we reproduce the unpolarized momentum distributions 
of Melnitchouk and Thomas jl^, , whereas the CS re- 
sults are inconsistent. The VNN vertex in Eq. ( |3.lD is 
also consistent with the one in Ref. . 

We also tested pNA helicity amplitudes in the vertex 
momentum (B), but the results disagree on the Jiilich 
expressions. However, if the momentum (A) is used, our 
results agree on them. It seems to us that the helicity 
amplitudes are shown for the choice (B) in pNN and 
for (A) in pNA. Therefore, as far as we investigated, 
we believe that the FS calculations are right also for the 
pNA process. 

In the following, we show the helicity dependent meson 
momentum distributions. Because the distributions with 
Xv = are irrelevant for calculating Af{y), so that they 
are not shown. The isospin factors are extracted out from 
the expressions. 



f\v f^,\ _ / ^'^1- ^VNB [y I '^iJ 



^ ' |i + -^iVTTT'- 1) j D^^^'^^ (,', ki) , 

oy y L yy m% J 

(Al) 



f-(fc^)„„x Jt,2 p2 ( I p \ 

f\v f^,\ _ / ^VNB yy ■ 

Jit, vnb \y) 



dy' y 



dy y' ^ 2ym^ ^^^^ ' ' 



n(fei)„ax 11.2 p2 / p \ 

J2L,VNB \y) 



(A2) 



dy' y 



m 



dy y' ^ y'ml ^^^^ 



(A3) 



/2T, VNB (y) 
dy' y 



VNB 



iy', kj 



k I m 



7 7^£Sr(v/^-i)^^'.^].(.',^i) 



dy y' ' 2y^y'mj^ 



(A4) 



Here, the partial derivative is given by 



dy' y_ 
dy y' 



-^2 ^ -1/2 

i + ^V(^i + "4) 

y mj^ 



(A5) 



The distributions Dy^"^ {y' , k^) are calculated for the 
prescription (A) as 

^VNN (2/' I ^i) 



gl {kl+y"m%) 



y'Hi~y'y 

+ gvfvy' [ki + y' {y'^ml - (1 - y'Wv)] 



f2 



/2 2 ; 2 I r '2 2 /I /\ 2 1 ^ 

y ■m^k^ + {y - [I - y )my \ 



(A6) 



2fc2 



y'Hi-y')'' 
J! 

4m', 



^VNN (2/'' ^i) ^ 

9v.fvy'{l-y') + -^{ki + y'^ml) 



(A7) 



^VNN iv'' ki) 

gvfv 



Xvk 



VKl_ 



2gly'^{l-y')mN 



y/3(l _y/)2 

+ -—1- [ki - y'2(i _ 2y')ml + (1 - y')^ml] 



TUN 
f? 



+ ^2y'mN {ki + y'^m% - (1 - y')ml} 



N 



(A8) 



In the VNA process, the distributions are calculated 
for the vertex momentum (A) as 



^VNA (y'' ki) 



■fvNA 



3y'3(i „ y')im%m'l 
ki+ki{3-4y'{l-y')}m\ 

+ ki[y'^2y'' + i2~y')^}m% 
+ 4j/'(l - y')^mNmAmy 

- 2y'2(l - y')^mNmlml + {1 - y')''m%m1. 



(A9) 



D^'mI {y, ki) - 



f2 

J VNA 



VNA\y^-±J 3y'3(i _ y/)2^^^2, 

; 4 2 I I 2 r /f '2 2 2 A 1/^ l\ 2 

k±m^ + k^{4:y rrif^m^ ~ Ay {I - y )niNmAmy 
+ (1 - y')^m*y} + 2,y'^mlm\ - &y'^{l - y')mNm\ml 



-I- 3(1 - y ) niAmy 



(AlO) 
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^vfvNA^-i- 



k\mi\] 



3t/''^(l — y')'^m\m'^ 

- 2y'ki {(1 - 2y')mNm\ + (1 - y')mAml} 

- (2 - 3y')y""^A^"^i - 22/"(l - y')m\ml 

+ 2y'(l - y')^m%mAmy - (1 - y')'^mNmy 



(All) 



In the same way, the distributions are obtained for the 
prescription (B) as 



D 



VNN - ^^3^^ _ ^^^2 

+ 9vfvy' {{l + y')kl + 2y''ml] 



9v [ki^+y m^) 



(A12) 



6y''^(\ ~ y'ym^my 
k\ {2mA + (2 - y')™^} + fci {V^A 

- 2(2 - 4t/' + ^/'2)TO^,m^ - 2(1 - y')m%mA 
+2(1 - + 2y'2mi - 2;'(1 - 2y'2)mwmi 

- 2j/'(l - y')mljm\ + 2y'{l - y')^m%m\ 



+ 2y'{l - y'fm%mA - y' [I ~ y'f 



'TV 



(A17) 



The fongitudinal distributions agree on the FS results 
in the hmit ^ oo except for a term in Eq. ( Al(^ ). 
The factor l/y'^(l — y'Y is written as l/?/'^(l — y'Y m 
Ref . [Q . It is possibly a misprint [psl . 



2fc2 



5^(1 -y')^ 



gvfvy'il -V') + ^ (^i + y"'"^) 



N 



(A13) 



Avfc I 



+ 



y/3(l_y/)2 

9vfv r,,2 /o Q„,/N„/™2 



ruN 



V{fci-(2-3y')y'm^} 



(A14) 



D 



L.+l 

VNA 



(y', fc! 



J VNA 



3y'3(l - y')^mj^m^ 
+ fci {(3 - 4y' + 4zy2)mi - 42;'(1 - y')^mNmA 
+ (1 - yTrnj,} + kl {y'^A - 4y' + iy'^)m\ 
- 2y'2(l - y'fmNm\ + 2y'{l - y'fmlm\ 
+ 4y'2(l - y'fm%mA ~ 2y'(l - y'fm%] + y 
+ 2y'3(l - y'^m^ml + y'^{l - y')^m%m 



..'4 6 
^A 



2y'^(l - y')'"^^< - 2^(1 - y'Tm^m 



+ y (1 - y ) "T- 



7V 



(A15) 



VNA 



^VNA (y ' k±) — T , 2 2 

3y''^(l — y')"^m^my 
+ fci {(3 + 2y')™A + 4y'TOjvTOA + (1 - 2y' + 2y'^)m%} 
+ kl {y'(6 + y')mi + lOy'^mjvmi 
- 2y'(3 - 4y' - y'^)m%ml - Ay'^l - y')m%mA 
+ y'2(l - y')^m%} + 3y'^m% + Gy'^niNmi 



3y'2(2 - 2y' - y'^)m%mi - 6y'^l - y')m%m 



+ 3y'Hl-y'fm%ml , (A16) 
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